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Abstract 



We determine the physical parameters (magnetic field and Doppler factor) of 
the homogeneous synchrotron self-Compton model allowed by the observed X~ray 
to gamma-ray spectra and variability of Markarian 501 during the 15-16 April 1997 
flaring activity. We find that magnetic fields between 0.07 G and 0.6 G and Doppler 
^ I factors between 12 and 36 could fit (depending on observed variability time scale) 

J^ ■ these observations. We take account of photon-photon pair production interactions 

of gamma-ray photons occurring both inside the emission region and during propa- 
k>>( , gation to Earth and find these to be extremely important in correctly determining 

V^ ' the allowed model parameters. Previous estimates of the allowed parameter space 

have neglected this effect. Future multi- wavelength campaigns during strong flaring 
activity, including observations from optical to TeV gamma-rays, should enable the 
physical parameters to be further constrained. 

keywords: galaxies: active - quasars: jets - radiation mechanisms: gamma rays - 
galaxies: individual: Mrk 501 



1 Introduction 

Four BL Lac objects have been detected in the TeV energy range: Mrk 421 (Punch et al. 
1992), Mrk 501 (Quinn et al. 1996), IE S2344+514 (Catanese et al. 1998), and PKS 2155- 
304 (Chadwick et al. 1998). In two of these, Mrk 421 and Mrk 501, the high level of 7-ray 
emission enabled the spectrum to be measured up to ~ 10 TeV (Zweerink et al. 1997, 



Aharonian et al. 1997, Samuelson et al. 1998, Djannati-Atai et al. 1998, Hayashida et 
al. 1998). In the case of Mrk 421, the spectrum can be adequately described by a single 
power law between 0.3 - 10 TeV. However the spectrum of Mrk 501 shows clear curvature 
over a similar energy range (Krennrich 1998). Recently, the spectrum of Mrk 501 has been 
measured up to 24 TeV by the HEGRA telescopes (Konopelko et al. 1999, Krawczynski 
et al. 1999). 

The 7-ray emission of these two BL Lacs shows very rapid variability. For Mrk 421, 
variability on a time scale as short as ~ 15 minutes has been reported (Gaidos et al. 
1996). In the case of Mrk 501, variability on a time scale of a few hours was observed 
during the 1997 high level of activity (Aharonian et al. 1998, Quinn et al. 1999 cited by 
Krennrich et al. 1998), and there is some evidence of variability on a time scale of 20 
minutes (Aharonian et al. 1998). 

The TeV 7-ray flares are simultaneous with the X-ray flares. In the case of Mrk 
501, during the 16 April 1997 flare the X-ray spectrum was observed by the Beppo-SAX 
observatory up to ~ 200 keV (Plan et al. 1998), and during the same high state OSSE 
observations made between 9 to 15 April 1997 (Catanese et al. 1997) show that the energy 
flux per log energy interval continues up to ~ 500 keV at roughly the same level. 

Gamma-ray emission from active galactic nuclei (AGN) is often interpreted in terms 
of the homogeneous "synchrotron self-Compton model" (SSC) in which the low energy 
emission (from radio to X-rays) is synchrotron radiation produced by electrons which 
also up-scatter these low energy photons into high energy 7-rays by inverse Compton 
scattering (ICS) (Macomb et al. 1995, Inoue & Takahara 1996, Bloom & Marscher 1996, 
Mastichiadis & Kirk 1997). In this model all the radiation comes from this same region in 
the jet. Such a picture can naturally explain synchronized variability at different photon 
energies. 

The inclusion of photon-photon pair production interactions of 7-rays with low energy 
radiation within the emission region has been used previously when constraining the 
physical parameters of blazars (e.g. Mattox et al. 1993, Dondi & Ghisellini 1995). This was 
also included in our previous paper (Bednarek and Protheroe 1997) where we constrained 
the allowed parameter space of the homogeneous SSC model, i.e. the magnetic fleld in the 
emission region and the Doppler factor, based on observations of the 1994 flaring activity 
in Mrk 421. More recently, a similar analysis has been performed for Mrk 501 by various 
authors (Kataoka et al. 1999, Tavecchio, Maraschi and Ghisellini 1998). However, in this 
recent work absorption on both the infrared background (IRB) and the internal radiation 
of the emission region has been neglected. In the present paper we show that inclusion of 
both these effects is vital in order to properly determine the allowed parameter space. 

In Section 2 we discuss the effect of photon-photon pair production interactions of 
7-rays with synchrotron radiation produced inside the emission region, and with IRB 
photons during propagation to Earth. In Section 3 we obtain the values of B and D 
allowed by the observed ratio of 7-ray to X-ray power, and in Section 4 we discuss 
additional constraints arising from the requirement that the radiative cooling time must 
be consistent with the observed variability time scale. In Section 5, we further narrow 



down the allowed model parameters by comparing the predicted TeV 7-ray spectrum with 
that observed, and finally we discuss the implications of our results. 

2 Absorption of Gamma- Rays by Photon Photon Pair 
Production 

For propagation of 7-rays through isotropic radiation the reciprocal of the mean interac- 
tion length for photon-photon pair production is given by 

^77(^7)" =—7^/ d^—T dssa{s), (1) 

where n{e) is the differential photon number density and a{s) is the total cross section 
for photon-photon pair production (Jauch and Rohrlich 1955) for a centre of momentum 
frame energy squared given by s = 2eE^{l — cos 6) where 6 is the angle between the direc- 
tions of the energetic photon and soft photon, and Smin = (2meC^)^, emin = (2meC^)^/4_E^, 
and Smax(£^, E^) = 4:8 E^. 

We shall now apply these formulae directly to obtain the mean interaction length in 
the IRB and, in a modified form, to obtain the optical depth in the synchrotron radiation 
produced in the emission region. 

2.1 Absorption of gamma- rays in the infrared background 

The spectra of 7-rays from extragalactic objects are modified by photon photon pair pro- 
duction interactions with the infrared background, and observations of the 7-ray spectra 
of Mrk 421 and Mrk 501 have been used to constrain the intensity of the infrared back- 
ground (Stanev and Franceschini 1998, Biller et al. 1998). Recently, Malkan and Stecker 
(1998) have modelled the infrared background by summing the contributions from differ- 
ent types of galaxy taking account of evolution. We use their upper and lower models 
which are consistent with recent COBE DIRBE data (Hauser et al. 1998) and COBE 
FIRAS data (Fixsen et al. 1998) which we use to extend the model of Malkan and Stecker 
down below 3 /xm. 

We obtain the mean interaction length, x.yy{E.y), in the IRB plus cosmic microwave 
background radiation (CMBR) and these are shown in Fig. |T](a). We also show separately 
the mean free path in the CMBR alone. Our results are in good agreement with those 
of Stecker and De Jager (1998) calculated for the same infrared background models (see 
Stecker and De Jager for references to earlier work). In Fig. |I|(b) we show the reduction 
factor, R{E^) = exp[—Tiii{E^)] = exp[—d/x^^{E^) for a source distance oi d = 202 Mpc 
obtained for Mrk 501 assuming Hq = 50 km s~^ Mpc~^ and a redshift of 2; = 0.033. 
Clearly, absorption in the IR background becomes very important above 400 GeV. Figure 
shows the 1997 April 15-16 HEGRA and CAT data together with the approximation 
used later in this paper for the high energy part of the spectral energy distribution (SED). 




The figure also shows this SED after correction for absorption in the infrared using the 
two reduction factors of Fig. |l](b). 

2.2 Absorption of gamma-rays in the blob radiation 

Absorption of 7-rays will also take place on the synchrotron photons produced inside the 
emission region. For the spectrum of these target photons we use a fit to the Beppo-SAX 
observations made during the April 15/16 flaring activity (Plan et al. 1998), together with 
an indication from OSSE Observations made during the high state in April 1997 (Catanese 
et al. 1997) that the spectrum continued to ^s.max = 500 keV with approximately the same 
energy flux per log energy interval. The differential photon flux in the optical to X-ray 
region observed from Mrk 501 during the 16 April 1997 flare (photons cm~^ s~^ GeV^^) 
is approximated by 

3 X 10-^£-^-^ e < Ss,! = 2.14 x lO-^GeV, 

2.5 X IQ-^e-^-^^ es,i <e< 6^,2 = 1-95 x lO^^GeV, , . 

1.86 X 10-6£-i-83 ^^^2 < ^ < ^^^g _ 2 X lO-^GeV, ^ ^ 

4.37 X 10-^£-2 ^^ .^ < ^ < ^^^^^^_ 

The photon density in the emission region depends on the dimensions of the source 
and the assumed Doppler factor. We assume that relativistic electrons are conflned inside 
a "blob" which moves along the jet with the Doppler factor D and has magnetic fleld B. 
In the homogeneous SSC model the radii of the emission regions of low energy photons 
(ri), X-ray photons (rx), and TeV 7-rays (r^) are the same. This region is constrained 
by the variability time scale observed, e.g. in TeV 7-rays, tvar, 

ri = r^ = rx ^ O.b^cDU^, (3) 

where .^ < 1, .^ = 1 being the usual assumption. The differential photon density in the 
blob frame of synchrotron photons is then given by 

where e = De' and e' are the photon energies in the observer's and the blob rest frames, 
and c is the velocity of light. 

The optical depth Tsyn{E') in the blob frame for 7-ray photons with blob-frame energy 
E' = Ej/D for e^ pair production inside the blob is given by 






rsyniE'^) = ^l de'^ '^'dssais). (5) 



We compute the optical depth by numerical integration of Eq. |^, and our results can be 
approximated very well in the energy range of interest, i.e. 0.1 TeV < E^ < 100) TeV, by 

T^y^iE,) ^ 3 X lO'D-'-'E^^-%^X-\ (6) 



for 1 < D < 100 and 10^ s < tvar < 10^ s. We note that the E^-^ energy dependence we 
obtain is consistent with equation (B7) of Svensson (1987) for a photon spectral index 
of 1.6. From Eq. ^ we find that D ^ 13 to avoid absorption at 1 TeV if the variabihty 
time is 20 minutes and applying .^ = 1. Clearly absorption by photon-photon interactions 
inside the blob can not be neglected. 

3 Constraint from Ratio of Gamma-Ray to X— ray 
Power 

The spectrum of Mrk 501 shows two clear bumps which, during the outburst stage, extend 
up to at least ~ 500 keV (Plan et al. 1998, Catanese et al. 1997), and up to at least ~ 10 
TeV (Aharonian et al. 1997, Krennrich et al. 1999, Djannati-Atai et al. 1998, Hayashida 
et al. 1998, see also Protheroe et al. 1998 for review), the highest energy of 7-rays observed 
being 24 TeV (Konopelko et al. 1999, Krawczynski et al. 1999). These multiwavelength 
observations of Mrk 501 allow us to define the ratio rj of the energy flux per log energy 
interval observed at a chosen 7-ray energy, E^ = 1 TeV, at which the emission is assumed 
to be due to Compton scattering, to that at a chosen X-ray energy, e = 2 keV, at which 
the emission is assumed to be due to synchrotron radiation, 

^ UE^dt^ '^^dedt ' UE'^dt' ^ '' de'dt' '' ^' 

where rtot(-E^) = Tsyn{E^) + ti^{E^) is the total optical depth for photons in the syn- 
chrotron radiation of the blob and the infrared-microwave background, and the primed 
quantities are measured in the blob frame. For the power at 7-ray energies we adopt the 
value reported by the CAT experiment at E"^ = 1 TeV (Djannati-Atai et al. 1998), and for 
the power at X-ray synchrotron energies we take the value at e = 2 keV (Plan et al. 1998). 
For these two energies rj ^ 3.6, and from this constraint we shall now obtain the required 
magnetic field as a function of the Doppler factor for various variability time-scales. 

The synchrotron spectrum at e' in the above formula (Eq. |^ can be approximated 
analytically by 

, dN , , dN , „ , „ 
e'^-T-de'^—dihy^{i), (8) 



de'dt dy 



where dN/d'y' is the electron spectrum in the blob rest frame. &syn(7') = kUsl' is the 
energy loss rate of electrons, where k = Acax/'i-, ctt is the Thomson cross section, and Ub 
is the magnetic field energy density. The characteristic energy of synchrotron photons is 
given by 

e' ^ 0.5£b7'^ (9) 



We use the observed low energy SED (Eq. |^) to estimate the blob-frame equilibrium 
electron spectrum 



-1.8 

"17 
dN 



di 



/— 1.8 f ^^ I 

ai7 7 < 7b,i> 

/-2.18 / ^1^1 

"27 7b,i <7 <7b,2, nn^ 

«37 7b,2 < y < 7b,3, 



I «47' ^ 7b,3 < 7' < 7b,max. 

where 7' is the Lorentz factor in the blob frame, 

7(,,n = Cley,^jDesfl\ (11) 

where n = 1, 2, 3, and ax = 7'b,i'^'^5 "2 = 1, 03 = 7'b,?, «4 = asl'h^s^, ^b = m^c^B/Bcr, 
Bcr = 4.414 X 10^^ G, rrie is the electron rest mass. 

The synchrotron spectrum emitted by electrons with power-law spectral index a, mul- 
tiplied by the square of the photon energy, is given by 

dN ,2 2akUBe" (2s'\-^''^'^^' 

e ^ K — , 12 



d£ dt c jj \^B/ 

and we use this formula for each power-law section of the electron spectrum. 

The ICS part of the Eq. |^ cannot be obtained analytically in the general case be- 
cause of the complicated form of the Klein- Nishina cross section, and so we compute this 
numerically using 

dN ^,2_ w2 r d^ r dNii',E'^) 



dE'df ^ ^A'. dy Je' . dt'de'dE. 

'Y 'mm ' mill 



-E'i = E't ^^ '^^^^de'di, (13) 
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where 7^, ^ E'^jm^^, e'^,^ = E'^/[4j'{j'-E'^/ 711^0% E'^ = E.JD, and dN{-f', E'^)/ dt'de'dE'^ 
is the ICS spectrum (see Eq. 2.48, in Blumenthal & Gould 1970) produced by electrons 
with Lorentz factor 7' which scatter synchrotron photons in the blob, and we use in this 
formula the soft photon spectrum given by Eq. ^ 

Having obtained formulae for the synchrotron power at 2 keV (Eq. 0) and the inverse 
Compton power at 1 TeV (Eq. [13|), we can substitute these, together with the optical 
depth reduction factors, into the equation for rj (Eq. |^). Setting r] = 3.6, we solve this 
equation numerically to obtain S as a function of D for various values of t^^r- The 
resulting allowed values of B versus D are plotted for ^ = 1 as the two thick full curves 
(for the low and high IRB models) in Figs. |^(a)-(c) for t^ar = 2.5 hours, 20 minutes, and 
2 minutes respectively. Fig. ^(d) is for ^ = 1/3 and tvar = 2.5 hours. In the next section 
we further constrain the allowed values of B and D by requiring the radiative time scales 
be consistent with the variability time scale. 



4 Constraints from Radiative Time Scales 

The simultaneous variability observed in TeV 7-rays and X-rays allow us to place a further 
constraint on the homogeneous SSC model. The observed decrease in the observed TeV 
7-ray and X-ray fluxes may only occur if the electrons have sufficient time to cool during 
the flare, 

t'cool < U.rD. (14) 

This condition is required if the model is truly homogeneous, i.e. the model we consider 
in the present paper which includes a homogeneous magnetic field and constant jet direc- 
tion and bulk Lorentz factor. In what follows, we consider the cooling time of electrons 
responsible for synchrotron radiation observed at 2 keV, i.e. 

7' = {2e/DsB)'^', (15) 

with e = 2 keV. The cooling time scale for synchrotron losses of electrons with Lorentz 
factor 7' is given by 

2 

/syn _ rrieC 
'^cool — 1 TT r ^ ' 

kUBj 

We next estimate the ICS cooling time. For the soft photon spectrum we adopt, some 
interactions with energetic electrons will be in the Thomson regime while others will be in 
the Klein-Nishina regime with relatively small energy loss. To calculate the cooling time, 
we neglect interactions in the Klein-Nishina regime, i.e. with photons above e'rp ^ rrieC^ /'-/', 
and use the simple Thomson energy loss formula 

,ics ^ rn^c' 



where 



f/rad(<4)= rn{e')e'de'. {U 

Jo 



Because the optical to X-ray energy fiux per log energy interval, Fe{E) = E'^F{E), 
increases with energy (i.e. the photon spectrum is fiatter than e"^) most of the inverse 
Compton energy fiux will occur near 

E^ ^ D-f'^e'^ = D-i'nieC^. (19) 

Thus, the ratio of the energy loss times for inverse Compton and synchrotron loss for 
electrons with energies such that they radiate synchrotron photons which are observed 
with energies of e = 2 keV is given by the ratio, p, of emitted energy fiux per log energy 



interval observed at e = 2 keV to the emitted energy flux per log energy interval observed 
at a 7-ray energy of 

E^ ^ {2eD/eBy^^meC^, (20) 

where we have used Eqs. |15] and |l^. The ratio is then 

Now, E^ = DE'^ depends on D and B, and rtot{E'^) depends on D, tvar, C, and B and the 
observed optical to X-ray SED, and so the ratio of emitted energy flux per log energy 
interval depends D, tvar, (,, B and the observed SED. 

The total blob-frame cooling time scale of electrons by both processes is given by 

/ 1 M~^ 

^ cool = I "Tics" + JiW^ ] ' \^^> 

\ '' cool '' cool / 

and this must be less than the Doppler factor multiplied by the observed variability time, 

111 

+ TTiyiT > 7^^- (23) 



J-/ICS ' j-/syn ^ p., 
'' cool '' cool -^f-v^. 

Using p{D, B, t^i,j., ^) we can rewrite this in two separate ways in terms of the ICS and 
synchrotron cooling times in order to investigate which process dominates the electron 
cooling for specific parameters of the blob, 

[1 + piD, B, t^ar, Ol^^var > <ool, (24) 

[1 + piD, B, U,„ 0'']DU.r > tZ^v (25) 

Each of the two equations above can be solved numerically for a given value of tvar to 
obtain a constraint on S as a function of D, and these two constraints have been added 
to Figs. |^(a)-(d). In each case, a separate curve is plotted corresponding to variability 
simultaneous with the TeV 7-rays occurring at 2 keV and 0.2 keV. In these figures the 
synchrotron cooling time scale of electrons is shorter than the observed variability time 
scale for parameters above the dot-dashed curves (given by Eq. ^5]), and the IC cooling 
time scale is shorter than the variability time scale for parameters to the left of the dotted 
curves (given by Eq. ^4|). Therefore, the parameter space allowed by the variability time 
scale condition (Eq. ^3]) lies above the dot-dashed curves to the left of the dotted curves. 
Before discussing the implications of these new constraints we shall consider one fur- 
ther constraint. This constraint arises from the condition that the maximum energy of 
electrons (determined by the maximum energy of synchrotron photons) must be higher 
than the maximum energy of 7-ray photons in the blob frame, i.e. 

7max"^eC^ > ^-^,max^"\ (26) 



with 7^ax — (2£^s,max-Bcr/''^eC^-B-D)^/^, and the maximum energy of synchrotron photon 
^s,max- This condition gives an upper hmit to the magnetic field in the blob 

5 < 4.4 X 10-2£,_^^^DE-Lx G, (27) 

where ^s.max is in keV and -E'^^max is in TeV. The low energy part of the spectrum from 
Mrk 501 extends up to at least 500 keV (Catanese et al. 1997) while the high energy 
part extends up to at least 24 TeV (Krawczynski et al. 1999, Konopelko et al. 1999). 
The constraint obtained using these values has been added to Figs. i|(a)-(d). In the 
next section we shall discuss the implications of the constraints obtained in this and the 
previous sections. 

5 Discussion and Conclusion 

We have now mapped out the allowed parameter space for ^ = 1 in Figs, ^(a)-(c) for 
tvar = 2.5 hours, 20 minutes, and 2 minutes respectively; Fig. ^d) shows results for 
^ = 1/3 and tvar = 2.5 hours. Allowed combinations of B and D lie on the thick solid 
curve corresponding to the infrared absorption model assumed (upper curve corresponds 
to lower IRB model). From the constraints on the inverse Compton and synchrotron 
cooling time scales, if the 2 keV X-ray fiux is observed to vary during the fiare with a 
similar time scale as the 7-rays, then allowed combinations of B and D must be to the 
left of the thick dotted curve labelled tij, and above the thick dot-dashed curve labelled 
ts,i. Similarly, if the 0.2 keV X-ray fiux is observed to vary during the fiare with a similar 
time scale as the 7-rays then allowed combinations of B and D must be to the left of the 
thin dotted curve labelled tiji and above the thin dot-dashed curve labelled ts,ii- Finally, 
the requirement that a 7-ray energy be less than the energy of the electron producing it 
constrains the allowed combinations of B and D to be below the thick dashed curve. 

Examining Figs. i|(a)-(c) we see that, as expected, as the variability time scale de- 
creases, the Doppler factor and magnetic field required both increase, and the allowed 
range of log-D and logS decrease somewhat. Comparing Fig. ^(d) for ^ = 1/3 with 
Fig. ^(a) for ^ = 1, both being for tvar = 2.5 hours, we see the effect of varying ^ is to 
increase the required magnetic field and to reduce the allowed range of logD. 

Not all of the parts of the thick solid curves within the "allowed" parameter space 
will actually give a viable SSC model as they may predict a high energy SED which may 
have a very different shape from that observed. To further delineate the allowed values 
of B and D we calculate the spectra of 7-rays emerging from the blob (after photon 
photon absorption in the blob synchrotron radiation) and propagating to Earth (through 
the infrared background). We have calculated the spectra for each of the points labelled 
A, B, C, or D, in Figs. ||(a)-(d) and normahzed these to the observed flux at 1 TeV. The 
results are plotted for the lower IRB model in Figs. ||(a)-(d) respectively. Fig. ^(e) shows 
the result of using the upper IRB model (points labelled A', B', C, and D' in Fig. |^a) 
instead of the lower one, and we see that a worse flt to the observed spectrum results. 



2.5 hr 


1 


lower 


0.1 


11.5 


20 mill 


1 


lower 


0.2 


20.4 


2 mill 


1 


lower 


0.6 


35.5 


2.5 hr 


1/3 


lower 


0.15 


15.8 


2.5 hr 


1 


upper 


0.07 


11.2 



Table 1: Parameters for the best fits to the Mrk 501 spectrum. 

Parameter Spectrum tvar ^ IRB model B (G) D 

Fig. 3(a) C Fig. 4(a) 

Fig. 3(b) B Fig. 4(b) 

Fig. 3(c) B Fig. 4(c) 

Fig. 3(d) C Fig. 4(d) 

Fig. 3(a) a Fig. 4(e) 



We now discuss which of these spectra are consistent with the observations. For the 
longest variability time scale considered, 2.5 hours (Figs. |a, |d and ^), the parameters 
corresponding to the points labelled A and B are ruled out because of their poor agreement 
with the observed spectrum, and point C gives the best fit, with point D being acceptable 
provided further absorption (e.g. due to photons from the accretion disk, other parts 
of the jet, or a dusty molecular torus, Protheroe and Biermann 1997) can steepen the 
spectrum. As we go to shorter variability time scales (Figs. |b and c) only points A are 
ruled out by the data. Table 1 summarizes the best-fitting model parameters for each 
case considered. 

As seen in Fig. ^(d), which is for t^ar = 2.5 hours and C, = 1/3, the constraint from the 
IC cooling time scale and the observed variability at 0.2 keV is very close to ruling out 
the best-fitting spectrum (point C), and definitely rules out the spectrum corresponding 
to point D. It has been reported by Buckley and McEnery (1999, cited by Catanese et al. 
1997 and Kataoka et al. 1999) that U band observations show variability simultaneous 
with the TeV 7-rays. Such U band variability would further constrain the allowed values 
of D as the corresponding constraint from the IC cooling time scale would be significantly 
to the left of the curve labelled tiji in Fig. ^(d), ruling out the best-fitting spectrum 
(point C). Inspection of Figs, ^(d) shows that spectra corresponding to points A and B 
in Fig. H(d) are ruled out by their disagreement with observation. Thus, the reported U 
band observations would rule out all models with ^ < 1/3 implying that the allowed blob 
radius is larger than ^cDt^ar/S. 

We would like to emphasize the importance of including photon-photon absorption in 
calculations when determining the allowed parameters of the SSC model. For example, 
in previous papers the energy at which the emitted energy flux per log energy interval 
maximizes in the 7-ray region is used when constructing constraints, as is the ratio of the 
peak luminosities. However, because of photon-photon absorption the energy at which 
the emitted energy flux per log energy interval maximizes can be significantly higher than 
the energy at which the observed energy flux per log energy interval maximizes, and this 
can lead to an incorrect determination of the allowed parameters. For example, we show 
in Fig. ^ that the maximum in the 7-ray spectrum shifts to lower energies by about an 
order of magnitude when photon-photon absorption is included; this has consequences for 



10 



the constraint obtained by Tavecchio et al. (1998) from a comparison of the location of 
the peaks in synchrotron and IC spectra. 

Finally, we point out the importance of simultaneous multi-wavelength observations 
which are vital to properly determine the physical parameters of the emission region. To 
be most useful, these observations should cover the full energy ranges of both the low 
energy (synchrotron) and high energy (Compton) parts of the SED. 
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Figure 1: (a) Mean interaction length for photon-photon pair production in the cosmic 
microwave background plus infrared background upper and lower models of Malkan and 
Stecker (solid curves). Dotted curve labelled 2.7 K shows mean interaction length in 
the cosmic microwave background, (b) Reduction factor, R{E^) = exp(— r) for a source 
distance of 202 Mpc. 
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Figure 2: Mrk 501 1997 April 15-16 HEGRA data (data binned in 1/10 decade intervals) 
and CAT data (data binned in 1/5 decade intervals). The error bars on the HEGRA data 
include systematic errors as given in Fig. lb of Krawczynski et al. (1999). The dotted 
curve shows the SED used in this paper, and the solid and dashed curves show the SED 
after correction for absorption using the lower and upper reduction factors of Fig. 1(b). 
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Figure 3: Constraints on the parameters (magnetic field B versus Doppler factor D) of 
the blob of assumed radius 0.5^ctvar-D for the 7-ray flare on 15-16 April 1997 in Mrk 
501 assuming different variability time scales of the flare: (a) t^ = 2.5 hr and .^ = 1, 
(b) t^ = 20 min and ^ = 1, (c) t^ = 2 min and .^ = 1, and (d) for tvar = 2.5hr and 
^ = 1/3. The full curves give the allowed values for B and D constrained by Eq. 7 for 
the case of absorption of 7-ray photons in the IRB using the lower model (upper curve) 
and the upper model (lower curve), the dot-dashed curves give the lower limit from the 
synchrotron cooling time scale (Eq. 25) and the dotted lines give the upper limit from the 
inverse Compton cooling time scale (Eq. 24) for electrons with energies which produce 
synchrotron photons with observed energies of 2 keV (thick lines marked by ts,i and ti^) 
and 0.2 keV (thin lines marked by ts,ii and tiji), and the dashed curve gives an upper 
limit from the maximum energies of synchrowon and inverse Compton photons (Eq. 27). 
The points labelled A, B, etc., mark the values of B and D for which we have calculated 
7-ray spectra and which we shall compare with observations in Fig. 4. 
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Figure 4: Gamma-ray spectra computed for the specific values of B and D using tlie lower 
IRB model are compared with the observations of Mrk 501 in the 15-16 April 1997 flare 
by the HEGRA (Krawczynski et al. 1999) and CAT (Djannati-Atai et al. 1998) telescopes: 
(a) t^ = 2.5 hr and .^ = 1, (b) t^ = 20 min and ^ = 1, (c) tj, = 2 min and .^ = 1, and (d) 
for tvar = 2.5hr and ^ = 1/3. The spectra shown in (e) are obtained using the higher IRB 
model, ty = 2.5hr and ^ = 1. Specific spectra corresponding to different values of B and 
D indicated in Figs. 4(a)-(d) by A or A', B or B', C or C', and D or D' are shown by the 
dotted curves, dashed curves, solid curves, and dot-dashed curves, respectively. 



16 



> 



\ -6 



o 




1 

\og{E^ / TeV) 



> 

tr' 



\ - 



O 




1 

\og{E^ / TeV) 



-5 


■ ■'■•: ' , . . . . 




CM 

1 




(e) ■ 


B 


\ •, 


■ 




^ ,j 


NJ 




M 






■--&.. 


■ 


> 








f 


1*JJ 


1 






■ 


(U 






^H^ 


*=■+ 


c 




" 


H 










^ 


-'^ 






. 


\ -6 


- 




H 


J-'\ 




- 


^^ 


\ 




\ 


\ 




, ^ 




\ 






vi 


\ 


^ 


\ 




c^ 






^' 


■ 


t^" 




\ 


] 


1 


. 


ao 






, > 


■ 


o 

1 — 1 


\ 




\>. 


■ 


-7 


\ 


.w. 


" 



1 

\og{E^ / TeV) 



Figure 4: (Continued) 
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Figure 5: The 7-ray spectra computed for the parameters B and D given by the point C 
in Fig. 3a. The dotted curve shows the 7-ray spectrum produced in the blob by relativistic 
electrons. The dashed curve shows the 7-ray spectrum modified by absorption in the blob 
synchrotron radiation, i.e. the spectrum emerging from the blob. The full curve shows 
the 7-ray spectrum after propagation through the IRB using the lower IRB model. 
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